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The problem of automatically protecting a quantum system against noise in a closed circuit is 
analyzed. A general scheme is developed built from two steps. At first, a distillation step is induced 
in which undesired components are removed to another degree of freedom of the system. Later a 
recovering step is employed which the system gains back its initial density. An Optimal- Control 
method is used to generate the distilling operator. The scheme is demonstrated by a simulation of 
a two level byte influenced by white noise. Undesired deviations from the target were shown to be 
. . . reduced by at least two orders of magnitude on average. The relations between the quantum version 

' of the classical Watt's Governor and the field of quantum information are also discussed. 
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O . I- INTRODUCTION 

I ■ Watt's Governor (WG), which was built in 1782, might be the very first machinery to solve automatically a control 
I problem. As an automatic control tool. Watt's governor aims to conserve some physical properties of a system 
subject to stochastic noise while maintaining its internal dynamics. Schematically, the WG can be viewed as a two 
T-H ' step process, measurement and correction. At the first step a measurement of the system is performed to check 
^ , whether the constraint has been violated. Next, if such a violation was found, a correction step takes place, and 
drives the system back to the allowed boundaries. 
I A quantum Governor is a natural requirement when the limit of nano machines is approached. Quantum computing 
I ' [1, 2, 3, 4] is another candidate for such a device. However, quantum mechanics imposes non trivial restrictions on 
I I the development of the quantum Governor. A measurement, which is a main feature of the classical WG intervenes 
iy~j . in the dynamics of quantum systems and it therefore should be avoided or reduced to a minimum. 

' A control scheme is traditionally categorized either as an open- loop or a closed-loop control [5, 6]. In the closed- loop 
control scheme one tries to extract information from a feedback from the quantum system in a way that allows the 
control of the system[7, 8]. Within this scheme, a controlled collapse of some of the wavefunction occurs and converts 
^ I ^ . some of the quantum variables of the system into classical parameters. One needs therefore to delicately balance the 
(— I amount of withdrawn information in order to conserve the quantum character of the system. Feedback control of 

I quantum systems has being extensively studied during recent years by several groups [7, 8, 9, 10]. 
^ An open-loop control scheme corrects the system without any measurement. In order to realize the QG a tool that 

distills quantum systems in some automatic fashion, has to be built. Distillation steps usually reduce the density of 
the system, and hence, in order to conserve the density of a controlled system, an extra-step to enrich the system and 
compensate for the losses is also required. 

In this paper we suggest and demonstrate a physical realization for an open-loop QG scheme. Our open-loop QG 
' is a two step routine (Cf. Fig. 1 ). The scheme starts from an initial state, for example a diatomic molecule in its 
ground electronic and vibrational state. This state is then disturbed by noise. At the first stage of action, an external 
field is applied to the system with the purpose of distilling the undesired components. The rejected components are 
moved to another degree of freedom of the system. In the present model to an excited electronic state. At the next 
stage, the freely evolving dynamics is set to enrich the initial state and return it back to the initial density. 

We will show that the distillation step can be performed by a particular unitary transformation. Recently, Optimal- 
Control-Theory (OCT) has been applied to find the field generating the unitary operation [11, 12, 13]. In the present 
paper the methods of refs. [11, 12] are generalized to obtain the control field for non-unitary operations under 
non-unitary time-evolution governed by the Liouville equation for open quantum systems. This control field is then 
employed to derive the distillation stage of the QG. This step is followed by a field-free propagation of the system 
leading to the desired reconstructed state. 

In this paper the QG model is presented and the tools required for its construction. In Section (II) a simple QG 
model is presented. Section (III) elaborates on the OCT mathematical considerations required to achieve the distilling 
transformation. Section (IV) presents simulations of the QG. A discussion and conclusion is presented in section (V). 
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FIG. 1: Schematic principles of the QG: (a) An initial state is propagated under the influence of noise, (b) The initial state after 
the distortion, (c) The undesired part of the state is removed to another degree of freedom by a controlled transformation, here 
to another electronic surface. The initial state is purified but its density is reduced, (d) A properly designed decay transfers 
population back to the initial state and reconstruct its density during a free evolution of the system. 



II. THE QG MODEL 



The Hamiltonian of the quantum Governor QG model is partitioned to: 

H = H'' + flnoise + Hg 



(1) 



where, H° is the free Hamiltonian of the system, H„oise is the stochastic noise, and Hg is the control part, governed 
by an external field. The automatic control scheme could be applied to complex quantum systems. The principles of 
such a control scheme will be demonstrated by a simplified model composed of a single two-level qubit with frequency 
u)i: 
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(2) 



with the two states denoted by 1 1 1)^ , [2)^ | . The two levels can represent for example two spin states or two vibrational 
levels in a diatomic molecule. The qubit is then infiuenced by an external noisy field: 



1 
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where /j. is the dipole moment and f{t) is a white noise function which obeys: 

(/(i)> = {f{t)f{t')) = (:6{t-t') 



(3) 



(4) 



The target of the control is to conserve an initial qubit state protecting it from the noise. We first will describe the 
route to build a QG for a particular target byte. This approach will then be extended to a general target byte. 



A. The Conservation of a byte in its ground state 



The state of the system is described by a density operator in the energy representation. The target and initial state 
are chosen as: 
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FIG. 2: Level scheme for a QG. (a) QG which conserves a byte on its ground state. A two level system with frequency lui is 
influenced by noise. An external field can be applied on the byte to couple it with another byte with frequency lu2- Spontaneous 
decay is allowed only as depicted on the figure, (b) QG which conserves a general byte. The lower (target) byte is composed 
of two degenerate levels, and so is the upper (control) byte. The external field coupled between the two bytes, while the 
spontaneous decay is allowed only as depicted. 



The propagation in time of p under the influence of the noise leads to an undesired population on the excited |2) 
state. To restore the state, the qubit is coupled to an auxiliary qubit with frequency U2, {|l)e ) |2)e}- The second 
qubit can be realized for example as two vibrational levels within the excited electronic state. The distillation step is 
achieved by applying the unitary swap transformation: 
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(6) 



The outcome of this transformation is that all the undesired population is transferred to the of the auxiliary e 
byte. This step cancels also the phase between the two states (see Fig. 2a). 

Note that the transformation is done under Liouvillc evolution which allow also non-unitary transformations to 
take place. A possible one step solution for the QG might be to leave the population in the desired state untouched 
while moving all the other population to the desired state, e.g, with the distiling operator 
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This will use heavily the uncontrolled, non-unitary components of the Liouvillian operators. Working at the weak 
dumping limit and demanding a fast transformation turn this idea to be unpractical. We are therefore forced to use 
a transformation which is close to unitary, and add another step for the completion of the task. 
The free evolution step is generated by the following Liouville-von Neumann equation: 

+ J^d{p) (8) 

Where Cd is a particular dissipative Liouville superoperator. 

Cd induces a selective decay between the vibrational states: — > and \2)^ — > |2)^. For example a decay 
induced by the transition state dipole between vibronic (vibrational+electronic) which has the above selection rules. 

After the distillation step, both the lower and upper bytes reside on the target state within the single byte subspace, 
i.e., with the correct relative population and phase. The decay step then restores the population from the upper byte 
reconstructing the density of the target byte. A cyclic application of this scheme on a quantum system purifies it 
reducing the influence of the noise to minimum. 

This process can also be considered as a perpetual preparation of the desired state. Note, however, that this desired 
state itself is never achieved directly in any of the steps. The governor acts by enforcing the conservation of the state 
by properly aligning the state on the directions of the target. 
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B. The Conservation of a byte in a general state 



We next study the preservation of a superposition state which has the form: 



|+) = a|l) +6|2L 



(9) 



with 
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The governor utilizes the unitary transformation: 
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This unitary operator transforms the basis set from the original basis , |2)^ , , |2)g| to a new basis set 

I 1 l+)e ' l~)e}- Now, after the noise influences the initial byte, the original byte can be distilled by the 
transformed form of the swap operator O^, Cf. Eq. (6): 
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So that under the operator O^: 
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It has been noticed [14, 15] that, non-degenerate qubits are very difficult to handle due to the relative; c;olicrent phase 
that develops under the free evolution. To avoid this problem we take the two bytes for the conservation of a general 
byte as a two couples of degenerate states. Initially all of the population is on the target byte. After the distilling 
transformation both the target and the auxiliary bytes arc, within the single byte subspace, on the desired state, i.e.: 
\+)g/g, with the right phase between the two states of each of the bytes. 

The transformation of Eq. (11) corrects the error in both population and phase. The principle of the correction 
is to move a relative error between two states within a single byte, to a relative error between two bytes. A decay 
step must then return back the population from the control byte to the target byte and annihilate the relative errors 
between the bytes. This task is achieved using the same assumption used previously with the allowed transitions 
|l)g \l)g and |2)g \2)g, while other possibilities are forbidden. The scheme for the general byte conservation is 
illustrated at the lower panel of Fig. 2. 
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III. OPTIMAL CONTROL THEORY FOR NON-UNITARY TRANSFORMATION UNDER 

NON-HAMILTONIAN DYNAMICS 



The quantum governor is achieved by the unitary transformations responsible for the distiUation. The next step 
is to find the external field that induces such a transformation. This task is achieved by an inversion process which 
starts from the unitary operator responsible for the distilling and determines the field. The present description follows 
the treatment of Ref. [16]. A target transformation for a N-level space is described by the N-by-N matrix representing 
the operator O. O is neither necessarily unitary, or orthogonal. Nevertheless, practically it cannot deviate too much 
from miitarity. Our target is to find the field that generates the transformation O at time t = T^^ansi independent of 
the initial state. 

The density operator is now decomposed into a sum of a complete basis set of operators in the Hilbert-Schmidt 
space. The complete set for a N-level system density operators contains 2^^ hermitian matrices of dimension N-by-N. 
A scalar product between two operators A and B in Hilbert-Schmidt space is defined as[4]: 

(a • b) =Tr {AtB} (13) 

The norm of an operator is therefore: |A| = Tr |a^a|. For density operators: < \p\ < 1, so that \p\ = 1 for a 

pure state, and \p\ = 1/N for the classical limit for N degenerate states. Note that under unitary dynamics the norm 
and the entropy of a density operator are conserved. This, however, is not true under dissipative conditions. 

For a complete base of operators the desired operation 6, changes each of the basis set operators to new 

target operator: 

^target ^ qqOq^ (14) 

Under unitary transformation the complete orthonormal set is transformed to another complete orthonormal set. 
This is not true for non-unitary transformations. A chosen functional for the optimization procedure should reflect 

deviations between the propagated operators and the target set of operators. The set of operators jQ^^"^^"*! does 

not conserve the initial norm. Therefore is better to deflne the functional by: 



F = V (g° • GJ^*"'*) (15) 



where the set is obtained by propagating the set of target operators ^Q^^^sety J^ackward in time. When the 

target is achieved F = 2^ . 

Two additional constraints arc imposed: 

1. The reverse time-evolution of the system is also governed by the Liouville-von Neumann equation. 

2. The total field energy has to be minimized. 

To meet these demands a modified functional is employed: 



I Xit)\efdt 



, (16) 



Bj are 2^ operator Lagrange multiplier, and \{t) is a time-dependent scalar Lagrange multiplier. An extremum for 
F is obtained by a variation of F with respect to 6Gj and the field. Following Ref. [16], the equations of motion for 
the reverse propagation of the Gj 's become: 
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with the initial conditions Gj{Ttrans) = G*'''^*'^*, and for the forward time-propagation of the Bj operators: 
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Parameter 


Value 


A 


0.06601 Hartree 




7.2449716268 x 10"* Hartree 


fiU!2 


5.3746313155 x 10"* Hartree 


'^trans 


l.OSps 




lO.Ops 




241ps 



TABLE I: Parameters used in the simulation: A is the electronic energy gap, u!j is the vibrational frequency for the j—th 
electronic state, Ttrans is the duration of the distillation transformation, V is the decay rate, and r/ree is the period of a 
correction cycle. 



and Bj{Q) = G°. Application of the Krotov's iterative method to obtain monotonic increase toward the objective at 
each iteration defines the field at each new iteration at time t: 

e"-(f) = e^'-(i) + XC{t) ^ ( [a, G,{t)] ■ B,{t)) (19) 

3 

C{t) is a time dependent envelope function, usually a Gaussian, and A is a strategy parameter. Large values of A 
will cause rapid changes of the field at each iteration. The difficulty in building quantum computing components is 
believed to grow exponentially with the number of bytes for systems which evolve under unitary dynamics. For non- 
unitary dynamics it is expected to be even worse. This expected result originate from the fact that under non-unitary 
dynamics there are more options for decoherence to occur. 



IV. SIMULATION AND RESULTS 
A. Protecting a byte on its ground state 

The application of the QG in protecting a target byte of the form of Eq. (5) is now demonstrated. The first step 

of the procedure is to calculate the field required for the generation of the distilling transformation Cf. Eq. (6). 
Two main timescales dominate the QG model: 1) Ttrans - the time duration for the field-derived transformation, and 
2) T/ree " the period of free propagation, where Tirans Tfree- The characterizing parameters of the four levels of the 
model are taken as the two couples of the two lowest vibrational states within the first two electronic states of the 
Na2 molecule. Table I summarizes the parameters used in the simulations. 

The operation are carried out in a 4-level space. Due to the decay of the upper byte, it can be assumed that when 
the controlled transformation is applied, all the population is at the lower byte. Hence, only a basis set for a single 
byte ,i.e., the lower byte, operations subspace within a four level system is needed. The time evolution generated by 
the Liouville Van Neumann equation was calculated by the Newton polynomial expansion method [17]. 

Fig. 3 displays the infidelity of the transformation, log(l — vs. the number of iterations, n is the number of 

basis states involve in the operation, here 2. The two inset panels of fig. 3 show the resulting fields and their Fourier 
transforms which generate an infidelity close to -12. The two peaks in the frequency domain are the result of the 
initial guess which initiated the optimization process: e{t) oc cos(Ai) where A is the vertical energy gap between the 
two electronic states. 

At the second step, the target byte — \\) ^ (Ij^ is propagated freely for a duration Tfree with the additional 
distilling transformation. The white noise function for the propagation Cf. Eq. (3) is modeled as: 

f{t) = Ns{t) (20) 

where N is the noise intensity and —0.5 < s (i) < 0.5 is generated from a uniform random distribution. The decay 
rates are taken to be T^^ = ^22- For comparison, three calculations were performed with the same noise parameters: 

• A reference propagation: the system is propagated with noise without any correction. The red line in Fig. 4 
represents the value of i? = 1 — (/5° • p{t)) versus time. A significant deviations from the initial byte develops, 
up to ii = 0.02%, due to the noise. The time-averaged deviation is about one third of the maximal deviation. 

• A partially controlled propagation: the distilling transformation is applied on the byte each Tfree, but no decay 
between the bytes was allowed. The black line in Fig. 4 displays the deviations of the density matrix from the 
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FIG. 3: Main frame: Infidelity of the transformation gate Od versus the number of the iterations. Upper and lower inset panels: 
time and frequency dependency of the field. 



target fp vs. time, for this case. According to the transformation O^; any remains of population on the upper 
byte will be transferred back to the lower byte and ruin the efficiency of the noise reduction. As expected, the 
growth of deviations from the target byte develops here almost at the same rate as in the previous case. The 
maximal deviation is 0.01% and the average is again one third of the maximal value. 

• Fully Controlled propagation: the system is propagated with both the distilling transformation and the decay 
period between the bytes. The blue line in Fig. 4 represents deviations of the density matrix from the target for 
the fully controlled propagation. It can clearly be seen that the full scheme works well. The maximal deviation 
under the fully controlled propagation is reduced to 7 x 10~^% and the average to approximately 6 x 10~^%. 
Figure 5 is a blowup of this line. The influence of the very frequent corrections is clearly visible. The system is 
more stable by two orders of magnitude on average under the QG scheme. 

Changes of some of the parameters of the system, e.g., taking r/^ee to be smaller, can stabilize the system even 
further. 



B. Conservation of a general byte and an investigation of the necessary features for a QG 

For a general byte case the normalized state was chosen as the target byte: 

1 



l+)„ = l+) = 



(|o)-*fo„|i)) 



(21) 



and bn randomly chosen to be equal to 0.23113851357429. The two bytes were chosen as two pairs of degenerate TLS. 
The optimal field for this transformation gate was converged to the same accuracy as at the previous section with 
approximately the same speed of convergence. 

The noise model influenced directly both the population and the phase between the states: 



Hr, 



-1 1 
1 1 



(22) 



[1] Several numerical tests were performed for the present case as well as for the case of unitary transformation under unitary time evolution 
of ref. [12]. Typically the difficulty to achieve the optimal field depended very weakly on the nature of the desired gate. 
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FIG. 4: The deviation of the byte vs. time for the uncontrolled (red), partially controlled (black), and fully controlled (blue) 
cases. The deviation is defined as 1 — (p° ■ p(f)). See the text for the parameters used for the simulation. 
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FIG. 5: A blowup of the deviation of the byte vs. time for the fully controlled system: A distilling transformation is apply 
every T/ree = 241ps. 

To gain more insight on the necessary features needed for a QG to work properly several simulation for mutated 
QGs are presented. Two conventions for the characterizing parameters are employed: the first one is identical to the 
previous one: 



(23) 



A scheme which conserves low values of R can be defined as a fully conserved scheme. The distilling transformation 
aligns the byte to the right direction, i.e., with the correct relative phase and population between the two states. A 
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FIG. 6: The deviation of a general byte vs. time for the cases of uncontrohed propagation (black), equal (red), different 
(green) and exchanged (blue) decay rates, and for the case where the decay channel is not directed back to the subspace of the 
controlled byte(yellow). The two panels display (a) R - the un-normalized and (b) J?„ - the normalized deviations, respectively. 
See in the text for a more detailed description of the various cases 
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TABLE 11: Deviations from a general target byte for five different test cases (see the text for details). The columns display the 
maximal and the average of the deviations R and Rn for the five cases. 

byte can deviated from the target byte in its norm but still conserve a high resemblance to it with respect to the 
correct phase and population. A normalized deviation is defined as: 



\p{t)\ 



(24) 



The upper and the lower panels of figure 6 represents the measures R and i?„ for five simulations with the same noise 
parameters. The results are concentrated in tablell. 
Five test cases are simulated: 



1. A reference propagation without any correction depicted by a black line with closed circles. R — R„ 



10 



fts-to - 




-20 - 



-30 



0.0 



5.0x10 

Time [a.u] 



1.0x10 



FIG. 7; Logarithm of the deviation of a general byte vs. time for the cases of uncontrolled propagation (black), equal decay rates 
(red), non-degenerate byte (green) and scrambled transformation (yellow). The two panels display (a) R - the un-normalized 
and (b) Rn - the normahzed deviation. See in the text more detailed description of the various cases 



2. A fully controlled scheme with equal decay rate. 
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IQps, is depicted by the red line with the closed 



triangles. The efficiency of the QG is the same as for previous case of the ground state byte. Both R and i?„ 
are conserved to values well below the un-corrected propagation. The system is more stable by a factor of 20 at 
the worst case and by more then a factor of 200 on average. 

3. This scheme is identical to (2), but with different decay rates, T^^ — -^^22 = lOps. Only selective decay is 
allowed. The results are depicted by the green line with open triangles. The scheme seems to work at the same 
efficiency as the previous case. The relative phase that develops between the two bytes during the propagation 
due to their different energies is destroyed by the decay so that only the inner phase between the states is 
conserved. 

4. In this case the decay channels were switched so that the decay channels are |l)g ]2)g and \2)^ 

The results (blue line with closed squares), show that under this scheme the noise is accumulated and the un- 
normalized deviation from the target is constantly growing. However, the accumulated error is still bearable, so 
that the mutated QG manages to transform the undesired part to the other byte and _R„ is still conserved. 

5. In this test case (yellow line with closed diamonds), it was assumed that the population of the upper byte decays 
by some drain channel to a bath outside of the system. Under this scheme the remaining byte is well conserved 
but the norm of the state is reduced significantly. 

Another four test-cases simulations are presented in figure 7. The numerical results are summarized in tablelll. 
The lines with black circles and red triangles present as on the previous demonstration, the uncontrolled and fully 
controlled QG. The case presented by the lines with the blue squares has a non-degenerate upper byte. The non 
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TABLE III: Deviations from a general target byte for four different test cases (see the text for details). The columns display 
the maximal and the average of the deviations R and R„ for the four cases. 

degenerate states of the byte develops a phase that destroys completely the correction scheme. It is interesting to 
note that just as in the previous example, a difference between the two decay rates does not cause a significant change 
in the efficiency of the QG. 

The last case that was checked (in the hne with yellow diamonds), is a scrambhng of the upper byte transformation, 
which makes a Hadamard transformation to the upper byte consecutive to the regular distilling transformation. The 
result of this mutation is a bearable accumulated error, which is moved constantly to the upper byte and leaves the 
lower byte close to the target state, but with lower norm. 

V. DISCUSSION AND CONCLUSIONS 

An integral and crucial part of quantum computing and information research is devoted to quantum error correction 
(QEC) [18, 19, 20]. The main question in QEC is the following: Suppose A is sending quantum information to a 
receiver B. An unavoidable influence of a noise may distort the quantum information with probability p < 1/2. 
How would B be able to reconstruct the data that was sent to A? The solution to this problem is usually given by 
redundancy. Before sending his quantum information, A must duplicate his data in several copies from which B would 
be able to withdraw the original data to a very high accuracy. 

The task of building a quantum governor is close to QEC, but is different in both motivation and strategy. A QG 
main goal is to reduce the influence of noise on the channel between A and B, that is, to reduce p. Moreover, the 
strategy to achieve the control uses mainly the system itself and does not create extra information. Considering the 
fact that the scaling of the difficulty of building quantum computers is believed to be exponentially in the size of the 
system, it seems that the task of protecting a single byte from decoherence might be more important for quantum 
computing then the ability to use QEC. 

The target entitled here as Quantum Governor can be stated in two versions: 

• The frail QG - one or more of the system expectation values is constrained, e.g., energy, angular momentum 
etc.. This constraint is quite similar to the original constraint imposed by James Watt on his classical governor. 
The similarities between such QG and error corrections arc seemingly minor. 

• The robust QG - the task of interest is the full conservation of the state of the system, and not only one of its 
observables. This demand is more difficult then the one imposed by the Watt governor but it brings the robust 
QG closer to the error correction field. 

Due to the fact that uncomutative operators cannot be measured simultaneously, it is well understood that feedback 
control might be applicative to several kinds of frail QG, but surely not for any robust QG. 

In this paper the QG problem was solved for the robust case for a model two level system. In this case the 
difference between the two versions is not large. The robust QG problem is also related to the problems of refocusing 
and (dynamical) decoupling [21, 22, 23, 24]. Both approaches aim to reduce influences of noise coupled to a quantum 
system. However, the strategies of the two schemes are totally different. The QG problem is a state oriented problem. 
It demands the conservation of known states, from the influence of noise of an unknown form. A treatment of the 
noise in terms of stochastic quantum equations (see for example[25, 26, 27]) is therefore unnecessary. The problems 
of refocusing and decoupling are noise oriented problems. They try to immune unknown states from a noise with a 
known form. Accordingly, the solution to the last couple of problems, e.g., Bang-Bang methods and its derivatives, 
uses the known structure of the noise in order to build the appropriate immune decoupling. 

To summarize, in this paper the fundamental demands for the task of building a Quantum Governor were developed. 
The basis was set for a scheme to achieve automatic control on simple quantum systems. The scheme was demonstrated 
through simulations on two level system. A reduction of the noise by more then two orders of magnitude was achieved. 
The necessary features of a working QG under the present scheme were examined. 
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The extension of the present scheme to more complicated system requires additional study. Several other schemes 
and methods to achieve QG, for example the use of the Quantum Zeno paradox could be purposed. The exploration 
of these possibilities is still under investigation. 
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